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We develop a consistent semiclassical theory of spin dynamics for an isotropic ferromagnet with 
a spin S — 1 taking into consideration both bilinear and biquadratic over spin operators exchange 
interaction. For such non-Heisenberg magnets, a peculiar class of spin oscillations and waves, for 
which the quantum spin expectation value m = (S) does not change it direction, but changes in 
length, is presented. Such "longitudinal" excitations do not exist in regular magnets, dynamics of 
which are described in terms of the Landau-Lifshitz equation or by means of the spin Heisenberg 
Hamiltonian. We demonstrate the presence of non-linear uniform oscillations and waves, as well as 
self-localized dynamical excitations (solitons) with finite energy. A possibility of excitation of such 
oscillations by ultrafast laser pulse is discussed. 
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I. INTRODUCTION 

Magnetically ordered materials (magnets) are known 
as essentially nonlinear systems Localized nonlinear 
excitations with finite energy, or solitons, play an impor- 
tant role in description of nonlinear dynamics, in partic- 
ular, spin dynamics for low-dimensional magnets, with 
different kind of magnetic order. To date, solitons in 
Heisenberg ferromagnets, whose dynamics are described 
by the Landau-Lifshitz equation for the constant-length 
magnetization vector, have been studied in details, see for 
review Refs. In terms of microscopic spin mod- 

els, this picture corresponds to the exchange Heisenberg 
Hamiltonian, with the isotropic bilinear spin interaction 
JS1S2. 7 For a spin of S > 1/2 the isotropic interaction is 
not limited by this term and can include higher invariants 
such as (S1S2)™ with n up to 2S. In particular, the gen- 
eral isotropic model with the spin 5 = 1 and the nearest 
neighbor interaction is described by the Hamiltonian 

61 = ~ E [^(SiS^ + JTtoS,-) 2 ]. (1) 

<i,j> 

Here, the constants J and K determine the spin-bilinear 
(Heisenberg) and spin-biquadratic exchange interactions 
between nearest neighbors . This model ([1]) has been 
actively studied for the last two decades both in view of 
description of usual crystalline magnets, s ee pio neering 
articles Refs. and review articles Refs. llOHllI and in 
application to low-dimensional magnets, see recent Refs. 

For the model (fTJ) the character of the ground state is 
more complicated, than for Heisenberg magnets. It is de- 
termined by the values of the parameters of bilinear and 
quadratic exchange, J and K. In addition to the ferro- 
magnetic phase, which is stable at J > K, J > 0, and 



the antifcrromagnetic phase, which is stable within the 
mean field approximation at J < K, J < 0, two so-called 
nematic phases (collinear and orthogonal, see?^), are re- 
alized for this model. For these nematic states, the quan- 
tum spin expectation value m = (S) equals zero, even 
at zero temperature. (Below for short we will refer to 
the vector m as magnetization) . The areas of existing of 
the nematic phases separate from both sides the domains 
of stability for the ferromagnetic and antifcrromagnetic 
phases. Interest to the model (jXJ) increases in view of 
investigation of multicomponent Bose-Einstein conden- 
sates of neutral atoms with nonzero spins3^. At two 
chosen values of J/K, namely, at J = K and J = 0, the 
model ([I]) has the symmetry SU(3), which is higher than 
the rotational symmetry inherent to SO(3) ~ SU(2), and 
in a one-dimensional case is exactly integrable. The lat- 
ter is interesting from the theoretical point of view. 

A possibility to change the magnetization m in length 
is an important peculiarity of the ferromagnetic phase in 
the model ([T]). It is worth noting, for the regular Landau- 
Lifshitz equation, frequently employed for description 
of spin dynamics, the magnetization length keeps con- 
stant. This property is associated with the fact that 
the Landau-Lifshitz equations naturally emerge within 
the approach of spin coherent states, or states of the Lie 
group SO (3) ~ SU(2). They are parameterized by a 
unit vector, the direction of the latter coincides with the 
quantum expectation values for the spin operator m, and 
all quantum expectation values of some products of spin 
components are expressed through corresponding prod- 
ucts of expectation values for components of spin opera- 
tors (dipolar variables), see for review Refs fl5lfl6l . 

For a Heisenberg ferromagnet with a purely bilinear ex- 
change (K = 0) the approach based on SO (3)— coherent 
states is exact, whereas for the model ((T|) one has to 
take into consideration quantum expectation values for 
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all irreducible operators, which include not only dipolar 
variables m, but also so-called quadrupolar variables, bi- 
linear on components of spin operator S. In principle, 
such variables can not be reduced to the m only; for ex- 
ample, (Sx—Sy) can be non-zero even the values (S x ) = 
and (S y ) = 0. In fact, Hamilton dynamics of the vari- 
able m = |m| takes place, and a variable, canonically 
conjugated to m, is quadrupolar variables of a struc- 
ture mentioned above. Such dynamics, which is appro- 
priately called longitudinal , in principle does not exist 
for a Heisenberg ferromagnet while considered within the 
framework of the Landau-Lifshitz equation or the Hamil- 
tonian (P) with K = 0. 

The fast change of the length of the magnetization m 
is of great interest now. Thermal quenching of magne- 
tization length, caused by ultrafast (femtosecond) laser 
pulse, is known for different ferromagnets nearly 10 
years i 17 ' 18 ' 19 Non-thermal laser control of magnetization 
is also realized ) 20 ' 21 see for review Ref. l22|. Principal 
possibility of dynamical (besides heating) quenching of 
m = |m| till the values m » and ultrafast dynamics 
of the variable m is of great scientific and technological 
interest. 

For the model |T]) ID solitons 2 ^ and 2D topological 
solitons in the collinear nematic phase ; 12 ' 25 and near the 
SU©— symmetrical point,— have been studied. How- 
ever, non-linear dynamics for other phases, even for the 
simplest ferromagnet phase, has not been studied yet. In 
this work we investigate 2D longitudinal non-linear spin 
oscillations and solitons in the ferromagnetic phase of 
a non-Heisenberg ferromagnet within a consistent semi- 
classical description of the model ([TJ) . In the Section 2 
the equations for a full set of spin quantum expectation 
values obtained within semiclassical approximation and 
describing effects of dynamical quantum spin reduction 
are discussed. Non-linear "longitudinal" spin oscillations 
in such the system are found there. Interaction of cor- 
responding non-dipolar degrees of freedom with electro- 
magnetic field is also discussed in this Section. 

In next Sections we have demonstrated the presence of 
specific longitudinal solitons, for which direction of the 
magnetization vector m remains constant, however, the 
magnetization changes in length. Such soliton solutions 
are obtained in the framework of the semiclassical equa- 
tions, in a continual approximation (Section 3) and by 
analyzing of a discrete problem for a simple square lat- 
tice (Section 4). The Section 5 contains conclusions and 
discussions of results obtained, as well as some overview 
of open problem. The discussion of the possibility of exci- 
tation of longitudinal spin oscillations by ultrashort laser 
pulse are also present in this concluding section. 



II. MODEL AND ELEMENTARY 
EXCITATIONS. 

To develop a semiclassical theory describing a magnet 
with a spin 5 = 1 with Hamiltonian (fTJ) and to make al- 



lowance for the spin reduction on a lattice site, we intro- 
duce generalized coherent states of SU(3) group param- 
eterized by a three-dimensional complex vector u + iv, 
see Refs. [ijH 

|u,v) = («J+*«i) I*j> . ( 2 ) 

j=x,y,z 

where \tj) are three Cartesian states for spin S = 1, u 
and v are real vectors. With account taken of the nor- 
malization requirements and an arbitrarity of the total 
phase, the vectors u and v satisfy the conditions 

u 2 + v 2 = l,uv = 0. (3) 

In terms of the variables u and v, all irreducible quantum 
expectation values for spin S — 1 states, including the 
magnetization vector (S) = m and quadrupolar variables 
(SiSj + SjSi), can be expressed by the simple relations 

m = 2(uxv),(5 i S' fc + 5 fc S'i> =2{5 lk -u l u k -v l v k ). (4) 

For the ferromagnetic ground state, which is stable for 
J > K, J > 0, the value of |m| = 1, while the state 
is degenerated in the direction of m. It means that in 
the ground state |u| = l/y/2 |v| = l/y/2. At |m| = 1 
rotation of these vectors in plane perpendicular m, does 
not change the state of a system. However, for any |m| < 
1, states, which differ from each other by direction of 
u and v in the plane, are physically distinguished due 
to anisotropy of quadrupolar variables. As we will see 
below, the angle of rotation of u and v plays a role of a 
generalized coordinate conjugated to the magnetization 
length m = m. 

Dynamics of the variables u and v for given spin on a 
point i in a lattice are determined by Lagrangian, 12 

L = -2h^2 v i ' {dui/dt) — IT{u,v}. (5) 

i 

where TT{u,v} is the system energy, which coincides 
with the quantum expectation value of the Hamiltonian 
([T]) calculated with the coherent states $Z§. For a lattice 
discrete model, an expression for the energy IT{u,v} is 
given in Ref. [l2|- Based on this Lagrangian, we can eas- 
ily analyze both the linear and nonlinear dynamics of 
the ferromagnet. In particular, using the explicit form of 
the energy IT{u, v} proposed in^ we can readily obtain 
the spectrum of linear elementary excitations (magnons). 
This spectrum contains two modes. The first mode does 
not depends on the biquadratic interaction constant K. 
Its dispersion relation has the same form as for usual 
Heisenberg ferromagnet, e(k) = 4 J (1 — C (k)), where 
C(k) = (l/2)[cos(fcza) + cos(fc y a)], k = p/H, p is the 
magnon momentum, a is the lattice constant (below we 
will limit ourselves with two-dimensional square lattice). 
In the long wave limit ka <C 1, k = |k| the usual 
parabolic dispersion law appears, e ~ J(ka) 2 . The sec- 
ond mode describes the oscillations of the modulus of 
magnetization, m = |m|, coupled with some quadrupolar 
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variables. It is natural to call them longitudinal magnons, 
see below Eqns. pi [TDjl. 

It is difficult to analyze nonlinear dynamics of the vari- 
ables u and v, since we have to operate with four inde- 
pendent nonlinear equations rather than with two equa- 
tions for angular variables, as in the case of the usual 
ferromagnet. However, it is possible to show that the 
full set of nonlinear equations for the u and v vectors 
has a partial planar solution for which the magnetiza- 
tion vector changes its length only, m = me 3 and vec- 
tors u and v rotate in the perpendicular plane (1, 2), 
where ei and e2 are unit vectors in the plane, per- 
pendicular to the magnetization vector m = me 3 , e^, 
i = 1,2,3 present an orthogonal set of unit vectors. Be- 
low, we will restrict ourselves to the analysis of such 
planar solutions. For this solution, only three quantum 
expectation values are non-trivial, namely, the magne- 
tization (T3 = in = 2(uiV2 — U2V1) and two quadrupo- 
lar variables, cri = (S 2 — S 2 ) = u\ — u\ + v\ — v\ and 
o"2 = (S1S2 + S2S1) = —2u\U2 — 2v\V2- One can eas- 
ily show that <j\ + a\ + a\ = 1; that is, the vector 
(7 = (7iei+ 0262+ 0363 is a unit vector. It is convenient to 
introduce the angular representation for this unit vector, 

03 = cos#, o\ = (S% — S 2 ) = sin6>cos(/>, 

0-2 = (S X S V + S y S x ) = sin 9 sin <f>. (6) 

The advantage of these variables is that they are un- 
ambiguously determined from a given physical state of 
the system. In contrast, the variables u and v contain 
the halved values of the angular variables 9, 4>, that re- 
flects the nature of vector u (or v) as a vector — director. 
Using the angular variables 0, </>, we can reduce the La- 
grangian ([5]) to the form 

z = ^( a0B e i -i).{^)-w{e,4 > }, (7) 

i 

where IU{0, 0} is the system energy, which depends on 
the discrete variables 9i, <pi. It is convenient to present 
the energy through the vector variable er, 

W = -\Ys [Kcr i a j + 2(J-K)a it3 a j , 3 }. (8) 

<ij> 

It is interesting to note, this Lagrangian formally coin- 
cides with that for a spin S = 1/2 uniaxial ferromagnet 
which is known as XXZ -model. For this model, the con- 
stant of isotropic exchange equals to K/2 and anisotropy 
of spin interaction is proportional to J — K. Thus, the 
general dynamics of SU(3)— coherent state for spin one 
magnet includes the particular class of solutions, which 
are described by classical model for a spin S = 1/2, that 
is quite unusual model for theory of magnetism. Limit 
cases correspond to the following simple physical mod- 
els: in the vicinity of a transition to the nematic phase, 
SU(3)— symmetrical point K — > J, an effective spin 
model ([5]) becomes isotropic, while at K/ J^Owe arrive 
at the Ising model. Naturally, anisotropy of the effective 



model is realized in the cr space and has not direct linkage 
to spacious rotations of spin operators. 

For the model ([7]) it is easy to obtain oscillations, which 
are a non-linear analogy of the above mentioned longitu- 
dinal magnons. For this excited states, at the lattice 
site \i the variable 9i — 9q = const not depends on time 
and 4>i = kli + tot. The frequency of such oscillations uj 
depends on the wave vector k and the amplitude 9q as 
following 

w(k) = (Acos9 /h){2(J-K) + K[l-C(k)}}. (9) 

Within the linear approximation, at 9q ~~ * 0, this fre- 
quency becomes the frequency of longitudinal magnons, 
previously obtained by Papanicolaou^ In the long-wave 
limit this spectrum becomes parabolic, and can be writ- 
ten as 

uj(k) = u cos9 (l + r 2 k 2 ), (10) 

where 

huj =8(J-K), r 2 a =Ka 2 /8(J-K). (11) 

Here the frequency luq is a gap of longitudinal magnons, 
ro determines a characteristic space scale. A st udy based 
on the Lighthill criterion, see for example Ref. [28|, shows 
that such uniform oscillations being excited in the system 
(for example, by ultrafast laser pulse, see below) are un- 
stable against self-focusing. As a result, essentially non- 
uniform states, like solitons, should appear. For their 
analysis it is easier to employ a continual approxima- 
tion, considering fa to be continuous functions of co- 
ordinates and time, 9i(t) — > 9(x,t),4> i (t) — > (p(x,t). For 
2D system, or for a thin enough film of magnet, which 
complies with standard geometry of experiment^ one 
can use 2D solutions, and present the Lagrangian of the 
problem as £ = j L{9, cf)} ■ d 2 x, where the density of La- 
grangian is 

L = {H/2a 2 )(cos 9 - l)(d<f>/dt) - w{9, <j>}, (12) 
the energy density w{9, <fi} is determined by the formula 

w{9, 4>} = (2/a 2 )(J - K) sin 2 9 + (l/4)Xsin 2 9 (V0) 2 + 
+ (1/4) [K + 2(J - K) sin 2 9] (V9) 2 . 

(13) 

As we will see below, solitons exist at < K < J, and 
we will limit ourselves to this region of the parameters. 

Let us discuss briefly an interaction of longitudinal de- 
grees of freedom with external fields, having in mind pri- 
marily a possibility of experimental excitation of such 
oscillations. First of all, a magnetic field H = He af- 
fects only the magnetization, for a planar solution it is 
m = (S) = e3 • cos6>. Therefore Zeeman interaction of 
<t with the magnetic field parallel to some direction e is 
described by the Hamiltonian H H — —H(e, 63)03. Actu- 
ally, the magnetic field directed in parallel with the mean 
spin m = (S) does not affect the system state, while for 
any other directions of H one can expect trivial change 
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of the orientation of unit vectors e, describing the pla- 
nar solution. It turns out, that ac- magnetic field is not 
effective for excitation of longitudinal oscillations. 

It appears, that excitation of longitudinal oscillations 
may be done through application of an electric field on 
the system. For simplicity, we will start with dc-electric 
field. Interaction of such electric field with a spin system 
of a magnet can be described phenomenologically on the 
basis of the following Hamiltonian 

jjiint) = lto*»)E. E (14) 

where £^ p is a spin-dependent part of dielectric per- 
mittivity and E = Ee is the electric field, e 2 = 1, 
see Ref. [29|. In principle, £^ can include all spin 
variables, describing the system state and allowed by 
symmetryj 38 i 39 In our case, the components of £y 
should include the contribution from quadrupolar vari- 
ables, e^EiEj = e Q ■ (SiSj + SjS^E.E,. The dis- 
cussion of the microscopic origin of such interaction, in 
particular, the value of the constant E®, is far from the 
scope of this article. 

A possible role of quadrupolar interactions can be 
demonstrated by simple example. Consider an electric 
field, perpendicular to the magnetization in the ground 
state m = me^. It can be easily seen, the influence of 
such field is equivalent to the action of some "effective 
field" on the variable <x of the same form as the Zee- 
man interaction of usual magnetic field with usual mag- 
netic moment, ij( int ) = — H^er. Here the effective field 
H Q is described by 



H Q = — {[(e 1 ,E) 2 -(e 2 ,E) 2 ]e 1 +2(e 1 ,E)(e 2 ,E)e 2 } = 
= — — (ei cos a + e 2 sin a), (15) 

where a is a doubled value of the angle between the vector 
of electric field E and the direction of unit vector ei in 
planar solution. Then a simple analogy between action of 
the usual magnetic field H on the magnetization m and 
action of the field on the vector a becomes obvious. 

For consideration of ac-field, for example, electric field 
of electromagnetic wave (light) it is enough to replace 
(l/8w)EiEj by (X/l6ir)EiEj , where Ei now is the com- 
plex amplitude of the time-dependent electric field. 29 
For linearly polarized light, the same expression (|15|) 
for the effective field appears. If light intensity is time- 
dependent, for example, for modulated laser beam or 
for ultrashort laser pulse, the effective field will be 
time-dependent. Being linearly coupled with variables 
oi , 02, it can excite the longitudinal spin oscillations 
found above. 



III. LONGITUDINAL SOLITONS IN A 
CONTINUAL MODEL OF A MAGNET. 

As Lagrangian does not depend directly on <f>, but on 
its derivatives, the model (|12ll3p has an integral of mo- 
tion, which determines the total spin projection on some 
axis. The same integral of motion is present for discrete 
model (|7I8[) . This integral of motion is suitable to be pre- 
sented via the number of spin deviations in the system 
N, for discrete model and its continuous counterpart it 
reads 

N = ^(1 -cosflj), oriV= [ (1 - cos 0)d 2 x/ a 2 . (16) 

In the framework of quantum mechanics, N possess inte- 
ger values and in line with Ref. we use this quantity for 
semiclassical quantization of solitons. The important in- 
tegral of motion is soliton energy E (J5J) which within the 
continuous approximation turns into E = j w{9, (f>}d 2 x. 

The integral of motion (fTB)) results in solutions with 
stationary dynamics, for which the vectors u and v, as 
well as in-plane components of the vector <x rotate in the 
plane (1,2) with some constant frequency. In this case, 
the variable 9 depends only on a distance from a certain 
point in the plane, which is considered as the center of a 
soliton. We limit ourselves to analysis of these solutions 
basing on ansatz of the form 

6 = 9(r),<p = wt, r = |x|, (17) 

Such solitons formally resemble so-called precessional 
solitons, which are known for the case of uniaxial Heisen- 
berg ferromagnets with precession of a unit vector of 
magnetization m, |m| = 1, around an easy axis (z-axis) 
with a constant frequency ujh and with an amplitude 
mj_ = sini9(r) depending on r, i?(r) — > at r — ► oo, see 
for review Ref. y. In spite of principal difference in phys- 
ical properties of these two types of solitons many formal 
features of them are similar. The latter allows us not to 
discuss some details. 

The function 9(r) is determined by a ordinary differ- 
ential equation, 

(6" + ^).(r 2 + ^sin 2 *) + 

lo ( 9' 2 a 2 \ 
+ — •sin0-sin0cos0 ( 1 — J = 0, (18) 

where the prime denotes the derivative over r, the charac- 
teristic size ro and the magnon gap frequency ujq are de- 
termined above (fTTjl . Far from a soliton the state should 
correspond to the ground state of the system, i.e. the 
condition 9{r) — > at r — > oo should be fulfilled. The 
condition 9'(0) = ensures the absence of singularity at 
r = 0. The equation (|18|l with such boundary condi- 
tions can be easily solved by the "shooting" method^ It 
has a discrete set of solutions 9 — 9 n (r) with n nodes at 
points r = r n ^ 0, n = 1, 2, .... Solitons with nodes are 
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unstable? 3 - and we will discuss the solution with n = 
and a monotonous decay of the function 9{r) only. 

Knowing the solution of 6{r) one can calculate integrals 
describing E and N , and represent the soliton energy 
as a function of the number of magnons bounded in the 
soliton, E = E(N). As already mentioned this procedure 
at N 3> 1 and within the continual approximation is 
equivalent to semiclassical quantization of solitons (for 
the discrete model some peculiarities occur, see Section 
5). It is convenient to use the fact that the equation l|18p 
can be formally obtained from variation of the functional 
L = E — HujN, 5L/59 = 0. The functional L coincides 
with Lagrangian ((Sj) calculated within the ansatz (fl~7|) . 
The condition SL/59 immediately leads to the relation 



hw = dE(N) /dN, 



(19) 



which coincides with that for precessional solitons. Eq. 
(fTn)l describes the quantum sense of the classical param- 
eter a; in the solution of the form (| 1 T|) : the value of ftw 
at N — > oo equals to a change of the soliton energy with 
a change of the number of bound magnons by one. 

Some limit characteristics of solitons can be obtained 
without an exact solution of (fT8|) . Using the phase plane 
method, it is easy to demonstrate that a soliton solution 
exists only at < lu < ujq and its characteristics depend 
on the parameter oj/luq. At (1 — ui/ojq) <C 1 the soliton 
amplitude is small and the function 9{r) takes the form 



9{r) = \J\ - ui/ujo ■ f (r/r )y/l - u/ojq 



(20) 



where f(x) is a universal function, localized in the area 
of Aa: ~ 1 with the value of the order of one in the origin. 
Further it is possible to demonstrate that for all K/ J at 
u> —>■ ujq soliton energy tends to a finite value E — > Eq = 
rj ■ AnK, where 77 ~ 0.93 is a numerical coefficient. In this 
limit the number of magnons is also finite, N — > A m j n = 
2tjN2, where N2 = 27r(ro/a) 2 is the characteristic number 
of magnons. 

These values are minimal for solitons in a magnet with 
a given parameter Kj J, and in the limit u> — » u>o the 
connection E — > fuo^N , which is typical for linear theory, 
appears. The similar property takes place for a preces- 
sional soliton with small amplitude^ with an essential 
difference, that for a precessional soliton the value of Eq 
is always of the order of the exchange integral J and can 
be compared with energy of a Belavin-Polyakov topologi- 
cal soliton Ebp = 47rJ, whereas for a longitudinal soliton 
for small Kj J the inequality Eq <C Ebp formally may be 
realized. In fact, for K <ti J the continual approximation 
is failed even at N ~ 1, and the minimal value of soliton 
energy E can not be smaller then 2.57 J, see the last para- 
graph of the next Section. On the other hand, the value 
of Eq found here within the continual approximation is 
valid for wide region of parameters like 0.7 J < K < J, 
where the value of N m [ n is larger then 4, see the next 
Section. 

Another limit case corresponds the condition < 10 <C 
ljq. To discuss it, we can mention that for u> = the 
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FIG. 1: The soliton energy E (in the unity of the exchange 
integral J ) dependence on the number of magnons N for the 
case K — 0.7J. The results of the continual approximation 
are represented by the solid line, the symbols are results of a 
numerical analysis of the discrete model, see the next Section. 
In the inset the detail behavior at small N is present; the 
arrow's end points out the end point of the dependence E(N), 
i.e. minimal value of soliton energy. The symbols beneath the 
arrow depict delocalized magnon states found by numerical 
simulation for a finite system. The dash line and doted line 
represent asymptotic behaviors at E — » huioN and square root 
asymptotic (|19[) for small and large N, respectively. 



2D soliton solution is absent, however, equations allow 
a ID "longitudinal" domain wall-like solution. For this 
solution, m(£) = m(£)e3, with m(£) — > 1 at £ — > 00 
and rn(£) — > — 1 at £ — > —00, £ is a coordinate along 
some direction in the magnet's plane. This wall has a 
characteristic width of the order of ro and energy a = 
a(K, J) per unit length. A qualitative analysis of (fT5|) 
gives that at < u> <C u>o a soliton contains a large 
enough circular region with a radius R ^> ro separated 
form the rest of the magnet by such a wall. Here again 
the situation is common to that for precessional solitons. 3 
Further it is easy to obtain a qualitative description of 
a soliton in this limit case. Apparently, that a uniform 
state with m(£) = — 1 has the same energy as at m(£) = 
1, and a finite region with m(£) = —1, i.e. 9 = n, does 
not contribute to energy, but affects the value of TV. In 
this case energy loss is connected only with the presence 
of a domain wall separating the inner region from the 
rest of the magnet. For the circular area energy loss is 
minimal at given N and one can see that N = 2wR 2 and 
E = 2ir<rR, where R is the radius of this area. Proceeding 
from that one can obtain a square root dependence of 
soliton energy on the number of bounded magnons N for 
a large soliton radius, which corresponds to the condition 
N> N 2 , 
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(21) 



Thus in the limit cases the dependence of soliton en- 
ergy on N, and also the parameters J, K is easily re- 
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FIG. 2: The same as in Fig.l, but for the case K — 0.5J. 
Please note essentially smaller values of iVmin and Eq 



constructed. In the intermediate frequency range, which 
corresponds to numbers of magnons N of the order of few 
N m in a thorough analysis of (fT8|) is needed. We carried 
numerical calculations for a set of values of K/ J for the 
region of interest < K < J. The analysis was done as 
follows: at the given K/J the equation (fig)) was solved 
numerically for a set of values of u> /uiq , which were chosen 
with different steps. Further having a solution the value 
of energy E and the number of magnons N were calcu- 
lated. Then the dependencies E(N) and u>(N) (the latter 
is important for analysis of stability of a soliton, which is 
stable in continual model at dio{N) / dN > only) were 
constructed. 

Let us briefly discuss the main characteristics of soli- 
tons found within continual approximation. The anal- 
ysis confirmed asymptotic dependencies derived above, 
see Figs. 1,2. In the whole region of parameters of the 
problem the function to(N) monotonously decreases with 
the N growth, i.e. the stability condition is fulfilled. 
Thus, in the framework of the continual approximation 
stable soliton solutions exist within the parameter region 
< K < J of Hamiltonian fl}. The soliton energies 
have a lower limit, Eq, which is smaller then the energy 
of familiar Belavin-Polyakov solitons. 



IV. THE DISCRETENESS EFFECTS FOR 
LONGITUDINAL SOLITONS. 

Strictly speaking, the continual approximation is valid 
only when the characteristic size of a soliton is essentially 
bigger than the interatomic distance, |V0| *C 1/a. This 
condition may be met for solitons with small amplitude, 
see (|20p. and also in the limit case J — K <C K, when the 
characteristic size ro is larger that the lattice constant 
a, ro 3> a. However, in contrast to precessional solitons 
in Heisenberg magnets with weak anisotropy, in which 
the characteristic length is tens and hundreds of the lat- 



tice constant, in our case the condition ro ^ a is much 
stricter. Even for enough small J — K = 0.1 • J, the value 
ro = 1.095 • a and only slightly exceeds the lattice con- 
stant a. In the region of the parameters K ~ J — K ~ J 
and in the especially interesting case K <C J, when the 
minimal energy of solitons is small, an applicability of 
this approach is not clear and one can expect essential 
discreteness effects. 

Let us consider the discrete model ([71 [5J for a square 
lattice. Analysis of discrete equations for the variables 
and <pi given for each lattice site demonstrates existence 
of a solution in the form of ifi — u>t and further it is 
possible to study only variables 0*. 

For analysis of solitons we employ the variation pro- 
cedure proposed and numerically realized in Ref. l32l . 
We will seek a conditional minimum of Hamiltonian, in 
fact, classical energy W(8i), with respect to variables 
9i, under the condition that the number of magnons 
N = 5Zj(l — cos#j) is fixed. While seeking a minimum 
one can find the precession frequency to from the equa- 
tion dW/d9i — hui ■ sin^i. It is worth noting, the sign 
of derivative du/dN in a discrete case is not important; 
a soliton is stable, if the found conditional extreme of 
energy is minimum. Analysis was done for an approxi- 
mately circular fragment cut from a square lattice sized 
24 x 24. We limit ourselves to such size as states we are 
interested in are essentially localized, and the influence of 
borders on them is negligible, while increase of a sample 
size require a significant increase of numerical calculation 
time. As one can expect at small values of (J — K)/K 
the behavior of the dependencies E(N) and ui(N) merely 
follows curves obtained within the continual approxima- 
tion, therefore we do not present them. As well, for a 
region of small N our analysis demonstrates that even in 
the case ro < a the results of the continual approxima- 
tion are quite close to numerical data, see. Fig. 1 and 
Fig. 2. 

It is interesting to note, that even for large N, when 
the characteristic size of inhomogeneity in a solution is 
of the order of r < a, these results qualitatively describe 
the dependence E(N) even at moderate values of K/ J, 
such as K/J = 0.7 and K/J = 0.5, for which r = 0.54a 
and ro = 0.354a, respectively. For K/J = 0.7 numerical 
data adhere closely continual curves, and a square root 
dependence with fitted value of domain wall energy a is 
working rather well. Even for the smaller value Kj J = 
0.5 only insignificant sign-alternating deviations from the 
square root dependence pTj) , which are almost invisible 
in Fig. 1 for K/J = 0.7, are observed on the numerical 
data on Fig. 2. 

However, apart from such characteristic of a soliton as 
E(N), which depends on integral values E and ./V only, 
effects of discreteness, nevertheless, are essential. This 
is apparently demonstrated in the dependence U)(N), see 
Fig. 3, 4 and especially it is clearly seen in analysis of 
the soliton structure, i.e. real distribution 9i in the lat- 
tice, see below. For dependencies ui(N) when N grows at 
the beginning, at small N regular deviation of numerical 
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FIG. 3: The dependence of the oscillation frequency u (in ljq 
units) on the number of magnons TV for a soliton in the case 
K = 0.7J. The results of the continual approximation are 
shown by the solid line, symbols are a result of the numerical 
analysis of the discrete model. The horizontal dash line points 
out theoretical value of magnon gap u)q. In the inset detail 
behavior at small TV is showed, the end of the arrow points 
out the end point of continual dependence TV = TV m j n . Sym- 
bols on the left side of the arrow correspond to delocalized 
magnon states, their insignificant deviations from the value 
of ujo are caused by boundary effects, naturally present for a 
finite system. 
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FIG. 4: The same as in Fig. 3, for the case K = 0.5J. The 
vertical arrows denote special values of TV = 18, 24 and 32, see 
the text. Please note that the dependencies correspond each 
other even at not small (u) — uq)/ujq < 0.25, that corresponds 
TV < 5. It is worth to mention that only one nonlocalized 
state with TV = 1 was revealed. 



data from the continual curve is observed. This feature 
can be explained by the fact that ui oc dE/dN, whereas 
the energy E(N) in the discrete model is lower than in 
continuum. However, even for Kj J — 0.7 it was observed 
noticeable deviations up and down from the smooth de- 
pendence typical for continual approximation, see. Fig. 




(a) TV =1 



(b) TV 




(c) TV =3 



(d) TV 



FIG. 5: Distribution of the discrete variables 9% and <f>i in 
a soliton at small -/V for the system with K/J — 0.5. Only 
a small part of fragment chosen for numerical calculation is 
presented. 9 and 4> f° r each spin are presented in the form 
of 3D vector cr = (sin 9 cos <f), sin 9 sin <f), cos 9). Values of 
5° < 9i < 90° are presented by arrows with open heads, 



< 



180 1 

:0 



by arrows with solid heads, the 



while for 90° < 

small values < 9i < 5 U are presented by open circles. The 
transition from non localized state at TV — 1 to localized states 
at TV > 1 is clearly seen. For TV = 3 the minimal value of 
|m| = |cos#| ~ 0.5 is observed in the soliton center. For 
TV = 8, the values of 9i nearly the soliton center are 158° and 
15°, and the tendency of transition to collinear structure is 
seen. 



3. For smaller value K/J = 0.5, this irregular behav- 
ior is much more essential, see Fig. 4. This complicated 
behavior is observed in that region of TV > 10 -C TVmin 
where the effective domain wall approximation and the 
square root asymptotic (|2Tj) should be applicable. There- 
fore as for precessional solitons in a Heisenberg magnet 
with strong single-ion anisotropy^ it can be naturally 
associated with characteristics of lattice pinning of a do- 
main wall. Let us discuss the structure of a soliton for 
K/J — 0.5, see Fig. 5 and 6. 

At small TV numerical analysis demonstrates almost 
radially-symmetric distribution of Bi with the scale of 
few lattice constant a, that is essentially larger than 
ro = 0.354a. It is interesting to note that even such 
a sensitive parameter as minimal value of TV is well re- 
produced by the continual calculation. According to our 
continuum calculation, for a magnet with K/J = 0.5 the 
value TV m i n = 1.6. The discrete analysis provides good 
localization of a soliton at TV = 2, 3 and much less local- 
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(a) TV =17 (b) TV =18 (c) TV =19 



FIG. 6: The distribution of the discrete variables 9i for a soli- 
ton for K/J = 0.5 at values of TV, close to "magic" value 
TV = 18. Values < 9i < 5° and 175° < 6 Z < 180°, corre- 
sponding to nearly maximal value of spin along the direction 
up and down, presented by open and solid circles, respectively, 
other notations and arrow scale here are the same as in Fig. 
|S] For collinear structure at TV = 18, as well as for TV = 24 
on Fig. [7] below, the mean spin deviation from the nominal 
value m = ±1 does not exceed 4 • 10 -5 . For adjacent values 
of TV, deviation from the values m = ±1 is not small, namely, 
value m ~ —0.75 in the corners of the soliton with TV = 17 
and m « 0.6 in the center of a wall segment at TV = 19. 



ized state at TV = 1, see Fig. 5. With further growth of TV 
up to N < 10 the distributions of 9i become much more 
sharp and tendency of formation of collinear states is ob- 
served. For the values of N > 10 the role of discreteness 
effects, primarily effects of domain wall pining, increases. 
This pinning may include either the dependence of the 
wall energy on its orientation related to lattice vectors or 
the dependence of the wall energy on its position related 
to distance to corresponding atomic lines in the lattice. 

A set of our numerical results may be explained con- 
sidering that, like in the case of uniaxial Heisenberg fer- 
romagnet, see Ref. HI, the most favorable position for a 
domain wall is to be placed between atomic lines like (0,1) 
or (1,0), but in contrast to Ref. [33, the domain wall is 
quite flexible and its bend from the line (0,1) to (1,0) does 
not costs too much energy. In principle, it corresponds 
to conclusions of Gochev, 31 who has demonstrated that 
for a 2D discrete classical spin model with anisotropy like 
(|5|) pinning effects are present for a domain wall paral- 
lel to the axis (0,1) and (1,0) and are absent for a wall 
parallel to atomic lines like (1,1). On the basis of these 
assumptions one can describe real distribution of the 0i 
amplitude and complicated behavior of u>(N) in a lattice 
soliton. 

In the case ro < a the favorable wall placed between 
adjacent lines like (0,1) or (1,0) is nearly collinear. There- 
fore the most favorable values of V are those for which 
a soliton contains a region where all spins have m = — 1. 
The region is separated from the rest of a magnet, where 
m = +1, by segments of such walls. This leads us to such 
spin configuration, for which this region would be of the 
form of a rectangular with the size l x xl y , where l x and l y 
are integers. Apparently, that the most preferable would 
be square areas with the spin number I and the magnon 
number N — 2l 2 . Following, 32 let us call these preferable 
states "magic" . Solitons with V = iV mag i c have the sym- 




(a) TV =23 (b) TV =24 (c) TV =25 



FIG. 7: The distribution of the discrete variables 9i in a soli- 
ton at TV, close to "half-magic" value TV = 24. Notations, 
arrow scale, and the value of K/J — 0.5 here are the same 
as in Fig. [6] The value of m » —0.75 in the corners of the 
soliton with TV = 23 and m » 0.66 in the center of a wall 
segment at TV = 25. 



metry axis C4 (hereunder we discuss the symmetry axis 
perpendicular to the lattice plane) , and spacious symme- 
try of a soliton coincides with the lattice symmetry. 

Along with magic numbers of magnons V mag i c = 2l 2 
"half-magic" states with a rectangular of spins in which 
l x , l y differ by one and TV = iVhalf-magic = 21(1 + 1) are 
also important. Their symmetry is lower then that for 
solitons with V = V mag i C and includes only the axis C2. 
For values like V = 21(1 + 2) no peculiarities were found, 
most likely they are always close to magic numbers, 1(1 + 
2) = (I + l) 2 — 1, and proximity effects to the latter are 
essential. 

Let us explain the character of lo(N) in a soliton with 
the value V, close to magic number. Configurations with 
N < V magic can be obtained by "corner smoothing" of 
a magic configuration due to wall band inwards "ideal 
square", see the example with V = 17 in Fig. 6. Since 
this does not require a lot of energy, the frequency in 
the region V < V mag i c is small and faintly depends on 
V, that corresponds well to our numerical calculation. 
If it is necessary to increase N up to V > V mag i c , then 
the situation is different: a domain wall should move 
into a unfavorable region of the lattice. An increase of 

V in the region V > V magic occurs only due to this 
not profitable wall sector and frequency values in this 
region arc quite large. It is important to note that for 
a purely collinear state (8 — 0. ir) the change of <j> does 
not change the system state, and the frequency value 
has no sense. Therefore at the transition of V through 

V = V mag ; c the uj(N) dependence experiences a jump. 
For this reason, two values obtained as frequency limit 
at TV = V magic + e and N = V magic — e are sketched in 
Fig. Q] In real calculation the value e = 0.01 was chosen. 

It is worth to discuss an important problem of soliton 
symmetry. The numerical analysis showed that in this 
case N > V mag j c the wall growth occurs only from one 
soliton side, see for example V = 19 in Fig. 6. Note 
the essential lowering of the soliton symmetry for such 
values of V, for which the soliton has no symmetry axis 
at all. The same regularities take place at the transition 
through half-magic value, however at N < Vh a if- mag i c 
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the soliton symmetry is lower, it contains an element of 
C2 rather than C4, see Fig. 7. Then, with increasing of 
N by values like 2-^3 the soliton symmetry is restored. 
Further when N reaches next chosen number (half-magic 
after magic or vice versa) this cycle reiterates, see Fig. 
4, where the positions of magic numbers N — 18, 32 and 
half-magic number N = 24 are depicted. 

In fact such tendency remains at extremely small val- 
ues of the ratio K/J up to K/J = 0.1 -r 0.3, when 
r ~ (0.12 -7- 0.23) • a < a. Again at N = 1, 2 
the soliton's size exceeds a, and the soliton can be de- 
scribed within the continuum approximation. Naturally, 
at small K/J solitons are more localized, for example, at 
K/J = 0.2 a purely collinear states appears at smaller 
N like minimal magic number N = 8. This quantitative 
difference leads to qualitatively new feature: at values of 
K < .KVcrit) — 0.36J even the states with N = 1 are lo- 
calized. These localized states resembles polarons, which 
are one-electron states localized due to interaction with 
non-linear media, see, for example^ These self-localized 
spin states with N = 1 can be called spin polarons. It is 
clear, the detail description of such states with small N 
should be based on exact quantum analysis, but we be- 
lieve, our semiclassical consideration gives at least quali- 
tative estimate for energies of such states. For such small 
K < if( cr it) <C J the maximal amplitude 9(0) is not 
small, the asymptotical solution (|20p is not valid more, 
and the energy of spin polaron state E sp = £A mm ) is 
smaller then the continual result 8(J — K). The value 
of spin polaron energy equals to 2.57 J at largest avail- 
able value of K = Kt CT n) ~ 0.36 J, and grows till 3.7 J at 

#«#(crit)- 



V. CONCLUSION REMARKS AND RESULT 
DISCUSSION. 

We considered spin dynamics in non-Heisenberg mag- 
nets with the spin S — 1 and biquadratic exchange taking 
into account. For such magnets, there are specific lon- 
gitudinal magnetic solitons, in the center of which the 
length of magnetization (S z ) is smaller than the nom- 
inal value (S z ) = 1, the values of (S x ) and (S y ) are 
equal to zero, but the oscillations of quadrupolar vari- 
ables (S x Sy + S y S x ) and (Sy — S x ) are presented. The en- 
ergy of such solitons is smaller then that of standard soli- 
tons described within the Landau-Lifshitz equation. In 
particular this energy is smaller then that for a "transver- 
sal" Belavin-Polyakov soliton, for which | (S) | = 1 and 
Ebp = 47rJ. Note that the Belavin-Polyakov soliton, as 
well as other transversal soliton states, are also presented 
in the model (fTJ) . 

Numerical analysis of the discrete lattice model demon- 
strated applicability of macroscopic approximation in 
vicinity of SU(3) point, where J — K -C J. For the 
rest region of parameters, this approximation is adequate 
for quantitative (at small N) or semi-quantitative, in 
general case, description of basic characteristics of soli- 



tons, for example, the dependence E(N). The analy- 
sis of the discrete model has also demonstrated a se- 
ries of qualitatively new effects, specifically important, 
when N is large enough and close to some chosen values 
of the magnon number, magic iV mag i C = 2l 2 and half- 
magic iVhaif- magic = 21(1 + 1), where I is an integer. For 
these chosen values a collinear state of solitons is real- 
ized, where the magnetization in all points has maximal 
values m = ±1. 

This result has been obtained within the semiclassi- 
cal approximation. We were not able to construct such 
states within exact quantum analysis of the model |T]). 
Let us discuss a role of quantum effects. An important 
question is whether or not the planar solution survives 
beyond semiclassical approximation. This problem can 
be discussed by taking into account the presence of usual 
ferromagnetic gapless "transversal" magnons within per- 
turbation theory. The existence of an exact semiclassical 
planar solution means that the equations for transversal 
variables, linearized over the planar soliton, have zero 
solution. In terms of magnons, this means that pro- 
cesses of single magnon radiation described by pertur- 
bation Hamiltonian H^t = Sfcl^p^p exp(iwt) + h.c] 
are absent (hereupon a p , aj, are creation and annihila- 
tion operators for such magnons with the linear momen- 
tum p and energy e(p) ~ J(ap/h) 2 ). In principle, there 
is a possibility for processes with radiation of several 
magnons, for example, two-magnon process describing 

(2) (2) t t 

by H lnt = Ei,2(*i,2 a i 8 2 ex pM) + ' i ' c '). where 1 = pi; 

three-magnon process, etc. Conservation laws of energy, 
momentum and total z— projection of spin, which can be 
written as E(N) - E(N - 2) = e(p) + e(-p), allow this 
process even at small E(N) - E(N — 2) ~ 2hu> < J, 
as the magnon dispersion law e(p) has no gap. One can 
expect the process of decay a soliton to magnons due to 
such radiations of magnons. As a result, the soliton will 
be characterized by finite lifetime. Such effects have been 
discussed early while going beyond the scope of semiclas- 
sical approximation for various topological solitons, 2D 
soliton with non-zero Pontryagin indexj22 and 3D soli- 
ton characterized by non-zero Hopf indexi^l However, for 
aforementioned examples these processes are slow and 
the lifetime of a soliton with large enough N is long, 
r - (h/J)N 5 and t ~ (h/J)N 5 / 3 in cases^ and^ re- 
spectively. Therefore one can expect that in our case 
with complete consideration of quantum effects solitons 
with large value of N will be enough long-lived excita- 
tions. Detailed discussion of this problem goes beyond of 
the scope of this work. 

For solitons in the discrete model one can point out one 
more interesting quantum effect absent at regular quan- 
tization of continual solutions with radial symmetry. For 
some special numbers of magnons N the lowering of soli- 
ton symmetry C4 inherent to the square lattice model ([T]) 
occurs down to C2 or even lower, see above Fig. 6 and 7. 
The presence of solitons with symmetry lower than lat- 
tice symmetry C4 means that in the classical case there 
are several (2 or 4) equivalent states, which differ from 
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each other by orientation in the lattice. In other words 
the classical state of the soliton is degenerated (two-fold 
or even four-fold) with regards to the soliton orienta- 
tion. In the quantum case there is a possibility for quan- 
tum tunneling (underbarrier transitions) between these 
states. For large N, the transition probability is low and 
can be calculated using instanton technique £>2£ As a re- 
sult one can expect the splitting of degenerated states, 
with creation of doublet or multiplet with four levels and 
lifting of the symmetry of the soliton to C4. We plan to 
return to detailed discussion of these effects in our future 
work. 

It is obvious that observation of effects of "longitudi- 
nal" spin dynamics is possible for materials with non- 
small biquadratic spin interaction. Yet Kittel demon- 
strated that such interaction appears due to interaction 
of spin system with lattice deformations.— For common 
reason, other mechanisms like electric multipole interac- 
tions and the Jahn- Teller effect equally resulting in bi- 
quadratic exchange (see, e.g., Ref. [33). There are a lot 
of such materials widely known, among them there are 
almost isotropic magnets, see review of NagaevJ^ 

In summary it is worth to discuss a possibility of exper- 
imental excitation of longitudinal nonlinear spin dynam- 
ics in the model ([1]) considered above. For a standard 
resonant method two problems come up. First, frequen- 
cies of these modes are rather high; second, magnetic 
field is coupled with dipole variables (magnetization) only 
and does not influence directly on quadrupolar variables. 
Both these problems can be solved by usag e of ultrashort 
intensive laser pulse, see, e.g., Refs. l40H4lT and for review 
Ref. H3- Usual value of a pulse duration r can be as short 
as 100 fs, and frequencies u> > 1/r, being considerably 
higher than frequencies of regular spin oscillations can be 



effectively excited. 

The possible role of different variables, dipolar and 
quadrupolar, can be demonstrated by a simple example. 
Consider a thick plane - parallel plate of a ferromag- 
net saturated along its normal (z axis). Let the light 
pulse propagates along the z axis, with the electric field 
parallel to the plate surface. The light interaction with 
dipolar degrees of freedom rtii = (Si) can be described as 
following. Due to inverse Faraday effect, a circularly po- 
larized light is equivalent to pulse magnetic field parallel 
to z axis i 40 i 41 Linearly polarized light produces two-fold 
anisotropy in the sample's planed Both scenario are in- 
effective for a sample, saturated along z— axis, and the ex- 
citation of usual transversal spin oscillations (magnons) 
is absent in this geometry. 

In contrast, the quadrupolar variables like (SiSj + 
SjSi) with i,j = x,y are coupled directly with linearly 
polarized light. The influence of such light pulse is equiv- 
alent to the direct action of some pulse of effective mag- 
netic field H^(t), see Eq. (fT5)) . on the variable <r of the 
form Jf( mt ) = — H^er. Being directed perpendicularly 
to the "ground state magnetization" er = e z , this pulse 
field effectively excite the oscillations of the x and y com- 
ponents of er, that is, the longitudinal spin oscillations 
considered in this article. The excitation of non-dipolar 
spin degrees of freedom by use of ultrafast optical pump- 
ing was recently observed for magnetic Mott insulator 
R 2 Cu0 4 42 
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